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ABSTRACT 

GRB 03 1203 was a very low apparent luminosity 7-ray burst (GRB). Coincidentally, it was also the first GRB 
with a dust-scattered X-ray halo. The observation of the halo allowed us to infer the presence of a large soft 
X-ray fluence in the total burst output. It has, however, also been claimed that GRB 031203 was intrinsically 
sub-energetic, representative of a class of spectrally hard, low-energy bursts quite different from other GRBs. 
A careful reanalysis of the available data, confirms our original finding that GRB 031203 had a very large 
soft X-ray component, the time of which can be constrained to within a few minutes after the burst, strongly 
suggesting that while GRB 03 1203 did indeed have a very low apparent luminosity, it was also very soft. 
Notions propagated in the literature regarding the uncertainties in the determination of the soft X-ray fluence 
from the halo data and on the available constraints from the hard X-ray data are addressed: the properties of the 
scattering dust along the line of sight (grain sizes, precise location and the geometry) are determined directly 
from the high quality X-ray data so that there is little uncertainty about the scatterer; constraints on the X-ray 
lightcurve from the Integral spacecraft at the time of the soft X-ray blast are not complete because of a slew in 
the spacecraft pointing shortly after the burst. Claims that GRB 031203 was intrinsically under-energetic and 
that it represents a deviation from the luminosity-peak energy relation do not appear to be substantiated by the 
data, regardless of whether the soft X-ray component is (arbitrarily) declared part of the prompt emission or 
the afterglow. We conclude that the difference between the soft and hard X-ray spectra from XMM-Newton and 
Integral indicate that a second soft pulse probably occurred in this burst as has been observed in other GRBs, 
notably GRB 050502B. 

Subject headings: gamma rays: bursts - X-rays: general - X-rays: ISM 



1. INTRODUCTION 

While 7-ray bursts (GRBs) are no longer as enigmatic as 
they were even a few years ago, the ability to use GRBs as a 
serious tool in cosmology and an understanding of their ba- 
sic mechanisms still elude us. Relations based on the energy 
release have the potential to resolve th ese difficulties. 

In particular, the Amati relation' (Amatietal. 2002|) be- 
tween the equivalent isotropic 7-ray total energy (£i so ) and 
the spectral peak energy £p ea k in G RBs, has been the fo- 
cus of considerable recent work (e.g. iBand & Preec e]l2005l 
iNakar & Piranll2005l iGhirlanda. Ghisellini. & Lazzatill2004 . 
Only a single burst apart from GRB 03 1203, has extended this 
relation to very low luminosit ies and peak e nergies (i.e. the 
low luminosity XRF 020903, Saka moto et all2003l) . 

It has also been suggested that the total energ y in 7-ray s 
from a GRB is nearly constant at ~ 10 51 erg (tFrail et all 2001). 
by correcting for the opening angle of the putative GRB jet. 
The determination of the opening angle is dependent on the 
time of the break in the lightcurve. This measure has proved 
difficult to use or understand because of 1) the difficulty 
in deciding the jet break time in lightcurves that are often 

1 Dark Cosmology Centre, Niels Bohr Institute, University of Copen- 
hagen, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark; darach, 
jens, pallja, kp @ astro.ku.dk 

- X-Ray Astronomy Group, Department of Physics and Astronomy, 
Leicester University, Leicester LEI 7RH, UK; sav2, rw, anl, pto, julo, jat, 
mgw @star.le. ac.uk 

3 Dept. of Experimental Physics, University College Dublin, Belfield, 
Dublin 4, Ireland: sfoley@bermuda.ucd.ie 

4 Laboratory for High Energy Astrophysics, Code 662, NASA 
Goddard Space Flight Center, Greenbelt, MD 20771, USA; 
jnr@milkyway.gsfc.nasa.gov 

5 Universities Space Research Association 



sparsely-sampled, contaminated by supernova features, and 
subject to fluctuations caused by density variations, and 2) the 
(few) cases where the total apparent energy release (equiva- 
lent isotropic) is well below this value. 

By any measure the apparent isotropic energy out- 
put i n GRB 031203 was extremely low (tWatson et alJ 
120041 hereafter W04), and for any opening angle of the 
jet, was significantly below the standard energy of ~ 
10 51 erg for GRBs inferred from jet op ening angles (W04). 
ISazonov. Lutovinov. & Sunvaevl (12004. hereafter SLS04) find 
an isotropic equivalent energy release of 4 ± 1 x 10 49 erg from 
the Integral 20-200 keV spectrum (an observed fluence of 
2.0 ±0.4 x 10 _6 ergcm- 2 ). Other bursts (e.g. XRF 020903) 
also have apparent energies below ~ 1 51 erg. 

It has been argued by SLS04 and bv lSoderberg et al J (120041 
hereafter S04) that GRB 031203 was representative of a new 
class of intrinsically sub-energetic bursts, possessing many of 
the characteristics of classical GRBs, but being a thousand 
times less powerful. This claim has far-reaching implications 
for GRBs. Ambitions to use GRBs as the most powerful dis- 
tance indicators in cosmology currently seem to lie mostly 
with the Epeak-E-y relation (similar to the Amati relation, but 
using the total co llimation-corrected 7-ray energy release, £ 7 , 
Ghirla r 1 d7etini200l . but whatever relation is used, a low- 
redshift calibration sample will be essential. If there is a dis- 
tinct population of under-energetic bursts, it will clearly need 
to be well-described and calibrated differently, especially if 
this type of burst dominates the low redshift sample. 

To suggest that GRB 03 1203 was intrinsically sub-energetic 
and a member of a new class of such bursts we must answer 
the question: was the total burst energy of GRB 031203 lower 
than expected compared to other GRBs? Such an apparently 
faint burst is expected to be soft according to the Amati re- 
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lation. Under the assumption that the emission detected by 
Integral comprised the entire burst, GRB 031203 was indeed 
much fainter than expected from this relation, since the In- 
tegral spectrum is hard. The high value of £p ea k adopted 
(> 190 keV), was based on the hard X-ray spectrum of the sin- 
gle pulse detected by the Integral satellite. But as we showed 
(W04 and IVauehan et al.ll2004l hereafter V04), the transient 
dust-scattered X-ray halo associated with the burst indicates 
that it was also very rich in soft X-rays, otherwise the halo 
observed by XMM-Newton could not have been so bright. 

The argument that GRB 031203 was a member (with 
GRB 980425) of a new, intrinsically under-energetic class of 
GRBs (SLS04; S04) hinges on the hardness of the burst. The 
fluence in the soft X-ray blast is critical to this discussion. 

The XMM-Newton data are therefore carefully reanalysed 
in this Letter. The dominant uncertainties in deriving the flu- 
ence are outlined in § ??. The full spectrum of GRB 031203 
and the consequences of analysing the complete dataset are 
presented and discussed in § ??. 

A cosmology with il m = 0.3, fl\ = 0.7 and Ho = 
75kms _1 Mpc" 1 is assumed throughout. Error ranges quoted 
are 90% confidence intervals, unless stated otherwise. 

2. METHOD AND UNCERTAINTIES 

Details of the XMM-Newton observations and the initial 
data analyses for GRB 031203 are outlined in V04 and W04. 
The luminosity of the soft X-ray blast, inferred from the dust- 
scattered halo observed by XMM-Newton, is key to the nature 
of GRB 03 1203. Here, we outline the procedure used to de- 
rive the fluence and analyse the major sources of uncertainty 
in this calculation. 

A complete model of the X-ray halo was used to find the 
best-fit parameters, including the rate of expansion, the width, 
the total fluence and the flux decay of the rings. The model 
produces a two dimensional distribution for the halo bright- 
ness with time and scattering angle for a given energy band. 

The fluence of the X-ray blast was inferred from the ob- 
served halo fluence divided by the scattering fraction. The 
differential scattering fraction as a function of scattering an- 
gle at a given energy is found by integrating scattering cross- 
sections over the dust grain size (a) distribution up to the max- 
imum grain size (a max ) and multiplying by the column density 
of dust. The uncertainty in the inferred blast fluence largely 
reflects the uncertainties in the scattering dust which is domi- 
nated by two things: 1) the size of the scattering dust column, 
and 2) the dust grain size distribution. 

2.1. The scattering dust column 

It was argued bv lProchaska et alJ J2004I) and later by SLS04 
that the fluence derived from the X-ray halo could have been 
overestimated by a factor of 4.4 in our previous work (V04). 
This was based on two incorrect assumptions. 

The first was that the individual rings observed in the halo 
were scattered by the total dust column along the line of sight 
(Ay = 3.6) 6 , whereas in fact the dust contributing to the rings 
is confined to relatively thin sheets of dust at well-defined dis- 
tances 7 1 395^3^ pc and 868^[gpc, see below). Dust that is not 
contained in these sheets cannot contribute to the scattered 

6 A higher Ay means a larger fraction of the X-rays are scattered which in 
turn would imply a smaller 'blast' fluence for a given observed fluence in the 
halo. 

7 The distance to the scatterer is known from fitting the halo's angular 
expansion with time; from geometrical arguments, D = 2ct/9 2 , where D is 
the distance to the scatterer, r is the delay between arrival times of directly 
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FIG. 1 . — Spectral energy distribution of the pulses detected using the dust- 
scattered X-ray halo and directly with Integral. The data from the X-ray halo 
are plotted as open circles with the uncertainty in the correction for the dust 
scattering plotted as a grey butterfly. The closed triangles represent data from 
the direct observation by Integral's IBIS instrument (SLS04), with the 90% 
limits to the best-fitting power-law (photon index, T = 1.63 ±0.06) plotted 
as dashed lines — the fluence at 1 keV derived from the halo cannot be made 
consistent with it. 



rings and since we use only the X-ray fluence in the rings 
themselves to derive the total fluence, other dust along the 
line of sight is irrelevant to the calculation of the burst flu- 
ence. Even using Ay = 3.6 as the extreme upper limit to the 
dust column contained in the sheets does not change our re- 
sults by more than a factor of 1.8. A t the same time it was 
also argued by Prochaska et al. ( 2004) that Ay ~ 1 of the total 
dust column actually belongs to the GRB host galaxy, which 
would leave only Ay ~ 2.6 as the upper limit to the dust col- 
umn available for the Galactic dust sheets. We find it unlikely 
that the entire dust column in this direction is contained in 
these two sheets. Based on the G alactic radial dust profiles 
(Neckel. Klare. & Sarcander 1980), the most likely value in 
the sheets is in fact Ay ~ 2.0 (V04). 

The second misapprehension was that the dust scattering 
fraction scales exponentially with Ay, whereas the depen- 
dence scales with the optical depth an d is therefore only lin- 
early related to the column density (Mau che & Gorensteinl 
1986). The factor 4.4 is the difference in optical extinction 
between Ay = 3.6 and Ay = 2.0, not the column density. This 
extinction relation is not correct for the X-ray scattering where 
the relationship is essentially linear. Since the maximum Ay 
has been argued to be ~ 2.6, the fluence in the burst could 
only have been overestimated by at most 2.6/2.0 which is 
~ 30%, though as noted above, this is unlikely. The effect of 
using Ay = 2.6 to derive the X-ray burst fluence is illustrated 
by dotted open circles and a lighter grey butterfly in Fig.^ 

2.2. The grain size distribution 

Because we possess time-resolved data for the evolution of 
the X-ray halo, the shape of the angular scattering response 
function for the dust (i.e. the way the scattered flux falls off 
with the scattering angle, Fig.|2) is strongly constrained. The 
largest grains always dictate this angular scattering response 

observed and scattered photons, and 9 is the observed angle, see V04 for 
more details. 
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FIG. 2. — Flux decay of the soft X-ray halo with time with best-fit model 
based on the dust-grain angular scattering profile. Over the period of the 
observation, the total flux in the scattered halo decreases as the scattering 
from increasingly large angles is observed. Increasing the grain sizes makes 
this decay rate faster because larger grains scatter more efficiently at smaller 
angles. The maximum of the grain size distribution (a ma x), dominates the 
scattering profile. We can therefore fit to the measured decay of the halo to 
find <3max- In this case, (3 m ax — 0.50 i 0.03 fiiu. 



function, allowing us to fit the differential scattering cross- 
sections to the observed flux in the halo as a function of scat- 
tering angle, with a max as a free parameter. This allows us to 
say that a max = 0.50 ± 0.03 /im along this line of sight. 

In our original analysis (V04) a single grain size was 
used. Here we assume a distribution proportional to a -3,5 , 
which gives a good representation of the optical extinc- 
tion curve and is simila r to that observed in X-ray scatter- 
ing for Galactic sources (|M athis 1 Rumpl., & Nordsieckl l 19771 
Mauc he & Gorensternll9 86: Prede hl & Schmitl 19951) 

Running the model with different values of the power law 
index of the grain size distribution, it is clear that values be- 
low -4.0 yield very large total scatteri ng fractions p er Ay, 
(> 12%), well above any observed value (Draine 2003). Even 
using the steep value of -4.0 implies a fluence only ~ 33% 
smaller than the our best estimate. 

The results from this analysis are consistent with our pre- 
vious estimate (V04) that used a single grain size, based on 
the dust scattering efficienc ies fou nd for Galactic X-ray halo 
sources (P redehl & Schmitlll 19951) . In other words, a similar 
Galactic source halo would have close to the same brightness 
for its central source as we infer for GRB 031203. 

The scattering efficiency is not very sensitive to variations 
in the details of the physical grain model. 

2.2.1. Dust scattering efficiency 

The dust model of Wein garjr^r&Draind (|2001 ) which has 
been used (e.g. Draine & Bond 2004) to convert optical ex- 
tinction (Ay) to X-ray scattering factor (r sca ), gives a scat- 
tering factor that is consistently 2-A times larger, over the 
0.7-3 keV range, than empirically determin ed from G alactic 
sources using ASCA, ROSAT and Cha ndra ( Woo et al JI 19941 
iPredehl & Sdmiit1[l99l iSmith et alJl200l lYao et alJl2003l) . 
Using this model, the fluence in the X-ray blast would be 
lower by a factor of ~ 3. However, as this model does not 
agree with the observational comparison of T sca and Ay (see 
Fig. 1 1 ofpDraine 20031), we have co ntinued to us e the empir- 
ically determined value from IPredehl & Schmittl (119951) . It is 
worth noting that the T sca -Ay relation is strongly dependent 



on the grain size distribution (see above), so that these are not 
independent sources of uncertainty. 

3. RESULTS 

Our reanalysis results in a 1 keV fluence density of 1320 ± 
260phcm" 2 keV- 1 (2.1 ±0.4 x lO^ergcm^keV" 1 ), a fac- 
tor of ten above an extrapolation of the Integral power-law 
spectrum to 1 keV (which has a 1 keV fluence density of 
110±20phcm- 2 keV" 1 and a photon index, T = 1.63 ±0.06, 
Fig. [0. The uncertainties quoted for the X-ray halo data in- 
clude the measurement error and the uncertainties related to 
the halo modelling. Given the size of these uncertainties and 
the fact that they are based on direct observation along this 
line of sight, we are forced to conclude that it is unlikely that 
the 1 keV fluence of the blast could have been substantially 
different. 

The analysis of the halo expansion was also improved by 
allowing the time of the X-ray blast to be a free fit parame- 
ter, Gaussian profiles were fit to the halo at different times to 
improve the radial size estimates, and the model fit was inte- 
grated over each time bin. We find results consistent with our 
previous work. The time of the blast was 600 ± 700 s after 
the beginning of the burst detected by Integral. The distances 
to the scatterers of 1395^0 pc and 868^[gpc are some of the 
most accurately known distances to any object beyond about 
50pc, with a total uncertainty of only < 2% at ~ 1 kpc. 

4. DISCUSSION 

The peculiar SED of the complete dataset, points to the fact 
that Integral and XMM-Newton observed different events in 
GRB 031203. A natural interpretation of these data is that 
there were two pulses in GRB 03 1203: one detected by In- 
tegral, with a hard spectrum peaking at or above ~ 190 keV, 
and a second pulse with a much softer spectrum, detected by 
XMM-Newton via its dust-scattered halo. 

It is expected that Integral's IBIS instrument should, in its 
lowest energy channels, have detected the harder X-rays as- 
sociated with such a powerful soft X-ray blast (Sazonov, priv. 
comm.). However, the lightcurve limits obtained by Integral 
are incomplete. Long (< 40 s) data gaps exist due to a ~ 100 s 
slew of the satellite. The slew occurred less than 300 s after 
the initial pulse. The IBIS data cannot therefore be used to 
place useful limits on the soft flux in the burst. (It may how- 
ever be used to place limits on the timing of the X-ray blast). 

Many bursts exhibit multiple pulses often accompanied 
by a strong softeni ng of the sp e ctrum, e.g. GRB 960720 
or GRB 97022 8 ( Front era et all l2000h or GRB 020410 
( iNicastro et all2004l) . The most striking case so far appeared 
during the preparation of this paper; the detec tion of a m assive 
soft X-ray flare in GRB 050502B (Burr ows etaHl 2005) start- 
ing ~ 500 s after the initial 7-ray pulse, and lasting ~ 500 s. 
The fluence in the soft X-ray flare was comparable to the first 
7-ray pulse, which had a hard spectrum. (Indeed, the photon 
spectral index, T = 1 .6, was very similar to that observed in the 
7-ray pulse of GRB 031203.) The consistency between the 
features observed in GRB 050502B and those inferred here 
from the X-ray halo of GRB 03 1203, reinforces the interpreta- 
tion that there were two very different pulses in GRB 03 1203. 

4. 1 . Afterglow or prompt emission ? 

The complete data show that although GRB 03 1203 was 
very faint (W04), there is no reason to suppose that it 
was anomalously so — it is more luminous than XRF 020903, 
for instance and probably of comparable luminosity with 



XRF 030723 llFvnbo et alJl200l iTominaga et alJl200l . The 
key issue is therefore whether it deviates significantly from 
the 'Amati relation', i.e. whether it was spectrally hard. The 
interpretation of the event — prompt emission, highly unusual 
afterglow, or reverse shock — while interesting speculation 
(SLS04), is not relevant to whether or not the burst was un- 
expectedly f aint. T he comparison is an observational one, i.e. 
Am ati et alJ §002) used the emission detected by the Bep- 
poSAX burst monitor and Wide Field Camera (WFC). To 
compare with these bursts in a meaningful way, we must 
use the same observational criteria and must include the soft 
X-ray blast in the calculation of the total luminosity, since 
its fluence or minimum flux would have been detected by 
the WFC dAmati et alJ l2002). The consideration of whether 
certain parts of the emission should or could be considered 
as afterglow or prompt emission is irrelevant for this com- 
parison, which is an observational one, based on the crite- 
ria for the sample selection. Based purely on the Integral 
data, GRB031203 appears to be one of only two signifi- 
cant outliers from this relation (the other being GRB 980425). 
However, when we include the XMM-Newton data, the X-ray 
(2-30keV) to 7 -ray (30-400keV) fluence ratio is S x /S 7 = 
1.8^9, which indicates that GRB 03 1203 was probably an X- 
ray flash, and certainly at least X-ray rich. This implies not 
only that the lower bound to the total X- and 7-ray fluence 
in the burst was roughly twice the 2 x 10~ 6 given by SLS04, 
but more importantly, that the peak energy of the total burst 



(if this is a well-defined concept in this case) was likely very 
low as originally concluded in W04. Taking this into ac- 
count, we conclude that there is no compelling evidence in 
the GRB energetics to suggest GRB 031203 was intrinsically 
under-energetic. 

In support of the argument that GRB 031203 was a cosmic 
analogue of GRB 980425, it was suggested that the shape of 
the prompt emission in was the same in both bursts (single 
pulse and FRED shape, S04). Since GRB 031203 could cer- 
tainly have possessed multiple pulses, this suggestion is not 
compelling. 

The luminosity of the X-ray afterglow at about one day 
is also very faint (9 x 10 42 ergs _1 at 10 hr, W04). The X- 
ray afterglow is, howeve r, still two orders of magnitude 
(Ramirez-Ruiz et al. 2005) brighter than predicted in the sub- 
energetic model proposed by S04 and both this and the low 
energy inferred from radio calorimetry (1.7 x 10 49 erg, S04) 
can be readily ex plained in a standard energy, off-axis model 
(Ramirez-Ruiz et al. 2005), which suggests an intrinsic peak 
energy for the burst of a few hundred keV, an order of magni- 
tude above the total observed XMM-Newton+Integral value. 
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